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Abstract 

Background: Spontaneous Bacillus anthracis mutants resistant to infection by phage APSOc (AP50^) exhibit a 
mucoid colony phenotype and secrete an extracellular matrix. 

Methods: Here we utilized a Roche/454-based whole genome sequencing approach to identify mutations that are 
candidates for conferring APSOc phage resistance, followed by genetic deletion and complementation studies to 
validate the whole genome sequence data and demonstrate that the implicated gene is necessary for APSOc phage 
infection. 

Results: Using whole genome sequence data, we mapped the relevant mutations in six APSO^ strains to csaB. 
Eleven additional spontaneous mutants, isolated in two different genetic backgrounds, were screened by PGR 
followed by Sanger sequencing of the csaB gene. In each spontaneous mutant, we found either a non-synonymous 
substitution, a nonsense mutation, or a frame-shift mutation caused by single nucleotide polymorphisms or a S 
base pair insertion in csaB. All together, S and 12 of the 17 spontaneous mutations are predicted to yield altered 
full length and truncated CsaB proteins respectively. As expected from these results, a targeted deletion or 
frame-shift mutations introduced into csaB in a different genetic background, in a strain not exposed to APSOc, 
resulted in a phage resistant phenotype. Also, substitution of a highly conserved histidine residue with an alanine 
residue (H270A) in CsaB resulted in phage resistance, suggesting that a functional CsaB is necessary for phage 
sensitivity. Conversely, introduction of the wild type allele of csaB in cis into the csaB deletion mutant by 
homologous recombination or supplying the wild type CsaB protein in trans from a plasmid restored phage 
sensitivity. The csaB mutants accumulated cell wall material and appeared to have a defective S-layer, whereas 
these phenotypes were reverted in the complemented strains. 

Conclusions: Taken together, these data suggest an essential role for csaB in APSOc phage infection, most likely in 
phage adsorption. (The whole genome sequences generated from this study have been submitted to GenBank 
under SRA project ID: SRA0236S9.1 and sample IDs: APSO Rl: SRS11367S.1, APSO R2: SRSl 13676.1, APSO R3: 
SRSl 13728.1, APSO R4: SRSl 13733.1, APSO R6: SRSl 13734.1, JB220 Parent: SRS1S0209.1, JB220 Mutant: SRSl S021 1 .1). 
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Background 

The evolutionary arms race between bacteria and their 
phages, giving rise to resistance and counter resistance 
mechanisms, shapes the bacterial and phage population 
size and composition as well as the evolutionary success 
of these two entities [1,2]. The major phage resistance 
mechanisms exhibited by bacteria and the strategies 
adopted by phages to counter these resistance mechan- 
isms have been extensively treated in two recent reviews 
[1,2]. Understanding phage resistance mechanisms is 
critical in developing phage therapy applications. 

Bacillus anthracis, a category A biothreat agent, is a 
spore-forming Gram-positive bacterium of the Bacillus 
cereus sensu lato group. It is a zoonotic soil bacterium 
that infects animals and occasionally humans, causing 
the disease anthrax. The October-November 2001 bio- 
terrorism attacks using mail laced with anthrax spores 
brought about renewed interest in research on B, 
anthracis. The prospect of biothreats using B. anthracis 
and the possibility of naturally emergent or deliberately 
created antibiotic resistant B, anthracis call for highly 
integrated and enhanced technological platforms for 
diagnosis and detection of this organism. This need is 
best illustrated in the case of a bacterial bioterror attack, 
where timely detection and intervention with counter- 
measures such as antibiotic therapy are paramount in 
preventing fatal consequences. Bacteriophages have been 
and still remain useful tools for bacterial species and 
strain differentiation [3-6], although evidence of success- 
ful application of phage therapy is still sparse in Western 
medicine [7]. 

Recently, the inherent binding specificity and lytic ac- 
tion of bacteriophage-encoded enzymes called lysins 
have been exploited for the rapid detection and killing of 
B, anthracis [8]. These authors demonstrated that the 
PlyG lysin isolated from the y phage of B, anthracis spe- 
cifically kills B, anthracis isolates and other members of 
the B, anthracis "cluster" of bacilli in vitro and in vivo. 
Both vegetative cells and germinating spores were shown 
to be susceptible. The lytic specificity of PlyG was also 
exploited as part of a rapid method for the identification 
of B, anthracis [8]. Thus, PlyG is a useful tool not only 
for the treatment but also for the detection of B, 
anthracis. 

The standard diagnostic tests for suspected B, anthracis 
as recommended by the Centers for Disease Control and 
Prevention (CDC) are as follows: presumptive identifica- 
tion to genus level {Bacillus family of organisms) 
requires Gram stain and colony identification. Pre- 
sumptive identification to species level {B, anthracis) 
requires tests for motility, lysis by y phage, capsule 
production and visualization, hemolysis, wet mount 
and malachite green staining for spores. Confirmatory 
identification of B, anthracis may include lysis by y 



phage, capsular staining, and direct fluorescent anti- 
body (DFA) testing on capsule antigen and cell wall 
polysaccharide. Thus, testing for y phage sensitivity is 
an integral part of B, anthracis identification [9]. y 
phage exhibits a fairly narrow host range but several 
B, cereus strains (e.g., ATCC 4342) have been shown 
to be sensitive to infection by this phage [8,10-12]. 
Several phages (CP51, CP54 and TP21) isolated from 
B, cereus and B, thuringiensis strains have been 
successfully used for transducing chromosomal markers 
and plasmids between B, anthracis strains [13-17]. 
However, their utility as B, anthracis diagnostic phages 
is limited because of their broad host range. 

Bacteriophage AP50 was isolated from soil using 
B, anthracis Sterne as the host [18]. Originally it was 
thought to be an RNA phage but was later shown to 
contain double stranded DNA and phospholipid [19]. 
AP50 was also shown to have a narrow host range; only 
one third of the 34 B, anthracis strains and none of the 
52 strains belonging to 6 different Bacillus species were 
susceptible to infection by AP50 [20]. Nine major struc- 
tural proteins were identified on SDS-PAGE gels. The 
molecular weight of the phage DNA was estimated to be 
9 X 10^ daltons [21]. Treatment with organic solvents 
such as chloroform (5%) and ether (25%) for 30 minutes 
inactivated the phage to a survival of about 1 x 10'^ 
[22]. 

In a recent study, we described the genetic 
characterization of AP50c, which is a derivative of AP50 
[23]. In contrast to the original report [20] on AP50, 
AP50c exhibited a much narrower host range, infecting 
111 of the 115 B, anthracis strains (97%) and none of 
the 100 B, cereus sensu lato strains. We also determined 
the genome sequence by 454 pyrosequencing and com- 
pared it to other Tectiviridae phages [23]. The genome 
size was determined to be 14,398 bp in length and 
shown to be highly similar to other Tectiviridae phages 
infecting Gram-positive bacteria in genetic organization 
and encoded proteins. We identified two mutations in 
the AP50 genome that were responsible for the control 
of lysogenic/lytic lifestyles and consequently the pheno- 
typic shift from turbid to clear plaque morphology. 
Based on these properties, we proposed that AP50c 
could be used as an additional diagnostic and thera- 
peutic tool for B, anthracis. Before embarking on phage 
therapy using AP50c, we sought to characterize the 
phage resistance mechanisms in B, anthracis in response 
to AP50c infection. In order to address this problem, we 
took an approach that began with whole genome scan- 
ning and concluded with genetic validation. Initially, we 
isolated a number of B, anthracis mutants that were 
spontaneously resistant to killing by AP50c and observed 
that these mutants produced an extracellular material 
that coats the cell surface [23]. We hypothesized that 
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this may be a bacterial adaptive resistance mechanism 
masking the bacterial receptor of APSOc. 

In the current study, we have mapped multiple bacter- 
ial mutations responsible for spontaneous APSOc resist- 
ance. Historically, mapping mutations responsible for 
specific phenotypes has involved transfer of the muta- 
tion with a closely linked selectable marker via one of 
the three conventional methods: transduction, conjuga- 
tion or transformation. Physical mapping has involved 
cloning and sequencing of the genetically mapped gene 
using traditional cloning vectors and methods. The re- 
cent development of rapid, cost-effective sequencing of 
bacterial genomes using second-generation sequencers 
or array-based comparative genome resequencing has 
accelerated the mapping of causal variations in bacterial 
as well as in other model system genomes [24-26]. We 
recently demonstrated the use of 454- generated draft 
sequences to map mutations conferring various pheno- 
types, and in that study a putative mutation associated 
with APSOc resistance was mapped to csaB [27], 
although the link between genotype and phenotype was 
not verified in that report. CsaB (for cell surface 
anchoring) protein is involved in pyruvylation of the 
peptidoglycan-associated polysaccharide, and this modi- 
fication is necessary for binding of the S -layer homology 
(SLH) domain. Thus, CsaB is critical for the non- 
covalent anchoring of SLH-domain containing proteins 
onto the cell surface [28]. In the current study, we ex- 
plore the relevance of the preliminary finding of a poten- 
tial linkage between the csaB gene and APSOc phage 
resistance by isolating and characterizing additional 
spontaneous mutants and, by using a recently described 
method [29], validating the WGS data by recreating csaB 
point and deletion mutations in a different genetic back- 
ground and showing phage APSOc resistance. Further- 
more, we restored sensitivity to APSOc infection in the 
mutants by complementing the defect in the AcsaB mu- 
tant in cis as well as in trans, thereby establishing an es- 
sential role of csaB in phage infectivity. 



Results and discussion 

Whole genome pyrosequencing and identification of 
putative variants in AP50'^ mutants 

We took a whole genome sequencing (WGS) approach to 
map the mutations and identify the gene(s) responsible 
for the phage resistant phenotype of the APSO^ mutants. 
We have shown previously that the draft sequences gen- 
erated by Roche/4S4 shotgun pyrosequencing can be used 
efficiently to map mutations responsible for various phe- 
notypes [27]. We generated whole genome sequence data 
of five spontaneous APSO^ mutant strains and compared 
them to that of the APSOc sensitive parent strain. We rea- 
soned that the causal variation of phage resistance might 
occur in the same gene (for example in the gene encoding 
the bacterial receptor of the phage) in all five mutants 
(S-Rl, S-R2, S-R3, S-R4, and S-R6), assuming that there 
is a single pathway to resistance. Hence, we searched for 
a common variation in all five mutants that is not present 
in the phage-sensitive parent strain. The phage-sensitive 
parent strain (Sterne 34F2) and the five mutants were 
sequenced by FLX or Titanium protocol. Some of the 
basic statistics of the genomic sequences are presented in 
Table 1. The estimated average coverage of the genomes 
ranged from 12X to 46X, and the number of high-quality 
variations (see Methods) ranged from 46-S2 (Table 1). 
The number of contigs was generally proportional to the 
Q39 scores; i.e., the lower the Q39 score the lesser the 
number of contigs. Among the high-confidence varia- 
tions, the most likely true positive variations (high-qual- 
ity) were selected for by applying the following criterion: 
variations that had a percentage concordance of more 
than 7S% variant base in the reads spanning any given 
position. Figure 1 shows the variations as a function of 
the percentage concordance of the variant reads. There 
were a total of 242 high-quality variations across all the 
five APSO^ mutants, and of these 2 IS showed 7S% or 
more concordance. Among these, 200 were identity-by- 
descent variations. The genomic locations of the IS 
unique variations along with their annotations are shown 



Table 1 Statistics of genome sequence of wild type and mutant B, anthracis strains 



Strain 


Runs 


Run type 


Length of assembled 
genome (bp) 


% < Q39 


Mean read 
length (bp) 


Number of 
contigs 


Depth of 
Coverage 


Number of putative 
HQ* variants 


34F2 


1 


FLX 


5,356,131 


0.06 


274 


101 


23 


48 


S-R1 


1 


FLX 


5,346,471 


0.37 


245 


227 


12 


47 


S-R2 


2 


FLX 


5,355,072 


0.02 


274 


60 


46 


50 


S-R3 




FLX 


5,355,679 


0.04 


283 


92 


28 


46 


SR-4 




FLX 


5,347,433 


0.06 


282 


102 


21 


47 


SR-6 




FLX 


5,349,716 


0.05 


277 


98 


26 


52 


JB220 




Titanium 


5,349,773 


0.39 


264 


111 


23 


43 


J-Rl 




Titanium 


5,348,852 


0.59 


269 


128 


22 


45 



*HQ, high confidence and high concordance (>75%). 
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Figure 1 Filtration of true positive variations from false positive variations by percent concordance and coverage. The graph displays 
260 high quality (HQ) variations in the genome sequences (34F2 and its derivatives) produced by GS Reference Mapper software as a function of 
percentage concordance. The graph indicates all the high-quality variations in the six genomic sequences obtained in this work. The variations from 
the six strains are indicated by different symbols. AP50 Rl, AP50 R2, AP50 R3, AP50 R4 and AP50 R6 are S-Rl, S-R2, S-R3, S-R4 and S-R6 respectively 



in Table 2. Among all the genes encoding variations, only Mutation in csaB is responsible for phage AP50c 

two were mutated in all five resistant strains and not in resistance 

the sensitive strain. These were BAS0840 (csaB), which Of the two mutated genes present in the AP50^ strains 

encodes a cell-surface anchoring protein [28] and identified by WGS, BAS0840 {csaB) was determined to 

BAS3946, which encodes a 2, 3-diketo-5-methylthiopen- be the most probable candidate gene for phage resist- 

tyl-1 -phosphate enolase. ance based on the following criteria: 1) In each of the 

Table 2 Location of true positive variations In genomes of mutant B. anthracis compared to reference genome 



Strain 


Start 
position 


Stop 
position 


Reference 
allele 


Sample 
allele 


Locus description 


Coverage 


Percent 
Concordance 


S-R1 


894,910 


894,910 




G 


CsaB protein 


10 


100 




3,894,626 


3,894,626 


G 


A 


2,3-diketo-5-methylthiopentyl-l -phosphate enolase 


8 


100 


S-R2 


210,679 


210,679 


T 


G 


intergenic 


29 


97 




895,097 


895,097 




G 


CsaB protein 


38 


82 




2,139,026 


2,139,026 


G 


T 


sensor histidine kinase 


34 


76 




3,894,626 


3,894,626 


G 


A 


2,3-diketo-5-methylthiopentyl-l -phosphate enolase 


24 


100 


S-R3 


895,097 


895,097 




G 


CsaB protein 


23 


100 




3,894,626 


3,894,626 


G 


A 


2,3-diketo-5-methylthiopentyl-l -phosphate enolase 


19 


100 


S-R4 


895,066 


895,066 


A 




CsaB protein 


18 


100 




3,894,626 


3,894,626 


G 


A 


2,3-diketo-5-methylthiopentyl-l -phosphate enolase 


19 


100 


S-R6 


160,743 


160,743 


A 




pXOl 


58 


86 




895,083 


895,083 




GC™ 


CsaB protein 


17 


88 




1,330,746 


1,330,746 


G 


CGGT 


intergenic 


4 


75 




3,894,626 


3,894,626 


G 


A 


2,3-diketo-5-methylthiopentyl-l -phosphate enolase 


14 


100 




4,303,652 


4,303,652 




AC 


solute-binding family 5 protein 


20 


100 


J-Rl 


68,981 


68,981 


A 




hypoxanthine-guanine phosphoribosyltransferase 


33 


97 




895,436 


895,436 


G 


T 


CsaB protein 


23 


96 
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Okb 2kb 4kb 6 kb 8 kb 

I BASOSag" ^ |BAS08 4Q> | BAS0841 | BAS0842 

csaA csaB sap eae 



csaB DNA sequence 



1 G S-Rl 

GTGCGGTTAGTTTTATCAGGATATTATGGTTTTTATAATGTTGGGG^ACGAAGCAATTTTGCAATCAATTATTAAAGCGCTACATGAAGAGGATCCT 

^ S-R8 

ACTCTTGAGCTTGTTGTACTTTCAAATGACCCCGACTATACAAGAAAAATGTACGGTGTAGAGGCAGTAAATCGCTGGGATATAAGAGCAATTTATA 

S-R4A S-R6 GCTTA S-R2,S-R3 G J-R2,J-R7 S-R9 T G J-R6 S-R7 G 

AAGAAATAAAGAGAAGTAATGGCTTA^ATTAGCGGAGGGGG^AAGCCTCCTTCAAGATAAAACGAGTATTAAAAGTATTTTGTACTATACAGGCATT 
T J-R4 S-R5 T 

ATGCGAATTGCTCGTTTTTTGAAGAAGCCGTATTATATTTATGCGCAAGGAATTGGTCCAATTACGAAAAGACAAAATCGTTTATTAGTGAAATGGC 

T J-R3 

AAGTATCAAAAGCAGAATATATATCAGTTCGTGATGAAGATTCATTTTTGTATTTAAAAGAAATTGGTATTAAGAAGGATATTGAACTAGTTCCAGA 

J-Rl T 

TCCAGTATTGGCATGTCAACCAGAAGGAATGAAGTCGGAGTGGTTGCAAAAACATTCGATTCAAGGGAAAGTAATTGCAGTTAGTGTGAGGTATTGG 
GACGCAAAAGAAGATTATATGAAGAAGCTAGCAGATACGCTGAAACAATTAAAGCGTGATGGATACCATATTTTATTCGTACCAATGCATGGACCAT 

C A693C-XinaI J-R8 A 

ttgatcaaaatgcatcacaooatattattaatttaatgooaoaagaagctcatatgcttccgtataagcttgatattcatgagaaaatttcaatttt 

j-r5 c gct 'h27 0a' 

otcggaatgttcgcttctaatcc;gtatgagacttcatgcgctcatactatctgccgttgctaatatccctatggtcggtatttcatatgatccaaag 
atcgattcatttttacaacaagtaaatcagccgattatcggaaacgtagatggtgattggacagctgaaactttgtataatgtcgcaacgaagcaac 
tagaacaaaaagaatatgtacaagaaacattggaacaaagagtagaagaattacgagaacaaatttcaacagcgtctcgatatatcataagtgactt 
gaattcaaaagagttcaaaaaaagaggaatgggatcttaa 

1104 

CsaB protein sequence 

Wild type s-ri s-r8 s-r4 s-r6 s-r2 j-r6 s-r? 

MRLVLSG YYGFYNVgTdEA I LQS 1 1 KALHEEDPTLELWLSNDpTdYTRKMYGVRAVNRVTOI RAI YKE iTkRSNGlTi SGGgTs LLQDKTS I KS I L Y 

J-R4I S-RSi J-R3I 

YTGIMfRIARFLKKPYYIYAQOIGPITKRQNRLLVKWyQVSKAEYISVRDTEDSFLYLKEIGIKKDIELVPDPVLACQPEOMKSEWLQKHSIQGKVI 

S J-Rl J-RSl J-R5 R A H270A 

AVSVRYWDAKEDYMKKLADTLKQLKRDOYHILFVPMHOPFDQNASRDIINLMGEEAHMLPYKLDIHEKlTsiLSECSLLIGMRLHALILSAVANIPM 
VGISYDPKIDSFLQQVNQPIIGNVDGDWTAETLYNVATKQLEQKEYVQETLEQRVEELREQISTASRYIISDLNSKEFKKRGMGS. 

367 

Mutants (Frame-shift and truncated CsaB proteins) 

S-Rl 

MRLVLSGYYGFYNVGGRSNFAINY . 
S-R8 

MRLVLSGYYGFYNVGDEAILQSIIKALHEEDPTLELWLSNDPTIQEKCTV. 
S-R4 

MRLVLSGYYGFYhT/GDEAILQSIIKALHEEDPTLELWLSNDPDYTRKMYGVEAVNRWDIRAIYKEIREVMA. 
S-R6 

MRLVLSGYYGFYNVGDEAILQSIIKALHEEDPTLELWLSNDPDYTRKMYGVEAVNRWDIRAIYKEIKRSNGLA. 
S-R2, S-R3, J-R2, J-R7 

MRLVLSGYYGFYNVGDEAILQSIIKALHEEDPTLELWLSNCPDYTRKMYGVEAVNRWDIRAIYKEIKRSNGLISGGGKPPSR. 
J-R4 

MRLVLSGYYGFYNVGDEAILQSIIKALHEEDPTLELWLSNDPDYTRKMYGVEAVNRWDIRAIYKEIKRSNGLISGGGSLLQDKTSIKSILYYTGIM. 
S-R5 

MRLVLSGYYGFYNVGDEAILQSIIKALHEEDPTLELWLSNDPDYTRKMYGVEAVNRWDIRAIYKEIKRSNGLISGGGSLLQDKTSIKSILYYTGIM 

RIARFLKKPYYIYAQGIGPITKRQNRLLVKW. 

J-R3 

MRLVLSGYYGFYNVGDEAILQSIIKALHEEDPTLELWLSNCPDYTRKMYGVEAVNRWDIRAIYKEIKRSNGLISGGGSLLQDKTSIKSILYYTGIM 

RIARFLKKPYYIYAQGIGPITKRQNRLLVKWQVSKAEYISVRD. 

J-R8 

MRLVLSGYYGFYNVGDEAILQSIIKALHEEDPTLELWLSNDPDYTRKMYGVEAVNRWDIRAIYKEIKRSNGLISGGGSLLQDKTSIKSILYYTGIM 
RIARFLKKPYYIYAQGIGPITKRQNRLLVKWQVSKAEYISVRDEDSFLYLKEIGIKKDIELVPDPVLACQPEGMKSEWLQKHSIQGKVIAVSVRYWD 
AKEDYMKKLADTLKQLKRDGYHILFVPMHGPFDQNASRDIINLMGEEAHMLPYKLDIHEKI . 



d 

AcsaB DNA and protein sequence 



Xmal 

GTGCGGCCCGGGGGATCTTAA 
M R P G G S . 



Figure 2 Genetic map of csa operon, DNA and protein sequences of csaB. a) Genetic map of the csaB region on B. anthracis chromosome. 
The blocl< arrows indicate the genes and the direction of transcription of the different genes, b) Nucleotide sequence of the wild type csaB gene 
and the location and nucleotide changes in the AP50^ mutants. Insertions and deletions are indicated by the red arrowheads and inverted 
arrowheads respectively. The various strains are indicated next to the mutated sites: derivatives of Sterne 34F2 and JB220 are indicated by the 
prefixes S and J respectively. The mutant loci labeled in black were discovered by WGS and those in green by Sanger sequencing of csaB PGR 
fragments. The altered bases are indicated in blue with the changed base above the wild type. A693G and H270A mutations were genetically 
engineered mutants. The sequences at the beginning and end of the gene marked in pink are the junctions the csaB deletion mutant, c) Amino 
acid sequence of wild type and mutant GsaB proteins. The sites of truncation of GsaB protein due to frame-shift and nonsense mutations are 
indicated by blue downward arrows. The amino acid changes in non-synonymous mutants are indicated in blue above the wild type residue. 
The amino acid sequences of the frame-shift and truncated proteins are shown. The residues in red are the additional residues added due to the 
frame-shift before truncation of the mutant proteins, d) DNA and protein sequence of the csaB gene in csaB deletion mutant. 
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five mutants (S-Rl, S-R2, S-R3, S-R4 and S-R6) for 
which WGS was performed, a frame-shift mutation (VG46 
in S-Rl, AA203 in S-R4, VGCTTA219 in S-R6, VG234 in 
S-R2 and S-R3) (Figure 2) was identified, which in all cases 
resulted in truncation of CsaB protein and possibly inacti- 
vation of its function. 2) Four additional spontaneous re- 
sistant mutants from the same genetic background (34F2) 
were screened by csaB PGR followed by Sanger sequencing 
of the csaB amplicons and were found to carry mutations 
in csaB as well (C385T, T274G, AC129, C254T). Similar to 
the findings in the first set of mutants, the C385T and 
AC129 mutations result in a truncated CsaB protein, 
whereas the T274G and C254T mutations result in a 
change from a noncharged (tyrosine) to a charged amino 
acid (aspartic acid) and a polar hydrophilic to a nonpolar 
hydrophobic amino acid, respectively. 3) A spontaneous 
AP50^ mutant generated in a different genetic background 
(JB220) was subjected to WGS and found not to carry the 
enolase (BAS3946) mutation, but it did carry a mutation 
in csaB (G573T), which results in a change from a posi- 
tively charged (arginine) to a noncharged (serine) residue. 
4) Seven additional spontaneous phage resistant mutants 
in the JB220 background were screened by the csaB PGR 
method described above and also found to carry mutations 
in csaB (G421T, C292T, G796G, T258G, C767A, and 
VG234, the latter of which was found in 2 different 



mutants). Interestingly, two of these mutations result in a 
noncharged residue being mutated to a charged residue, 
and the others result in a truncated csaB, These results are 
summarized in Figure 2 and Table 3. 5) Targeted introduc- 
tion of mutations in csaB in a strain (7702) not previously 
exposed to APSOc resulted in phage resistance (see below 
and Table 4). 6) Finally, BAS3946 was deemed unlikely to 
be involved in phage resistance because the mutation 
causes a non-synonymous change at G1024A leading to 
G342S in all five mutants (S-Rl, S-R2, S-R3, S-R4 and 
S-R6) and thus does not result in the loss of the protein. 
PGR amplification and Sanger sequencing of BAS3946 
from all the AP50^ mutants showed that the mutation is 
present only in seven out of nine mutants arising from 
34F2 background and the parent and not in any of the 
eight mutants of JB220 background. Additionally, a tar- 
geted deletion of BAS3946 or the G1024A point mutation 
in BAS3946 in a different genetic background (strain 7702) 
in a strain not exposed to APSOc had no effect on APSOc 
phage sensitivity (Figure 3). 

Reconstructed csaB point mutants and in-frame deletion 
of csaB result in mucoid colony morphology and phage 
resistance 

In order to unequivocally link the phage resistance 
phenotype to csaB and to rule out the involvement of 



Table 3 Results of Sanger verification of the various spontaneous APSO*^ mutants 



Sanger 
run ID 


Strain 


Sanger data 


Mutation 


454 data 


Implication 


BAS3946 


csaB 1 


34F2 parent 


WT 


None 


WT 




G1024A 


csaB 2 


S-Rl 


V G between positions 46 and 47 


Frame-sliift 


Same 


Truncation of CsaB 


G1024A 


csaB 3 


S-R2 


V G between positions 233 and 234 


Frame-sliift 


Same 


Truncation of CsaB 


G1024A 


csaB 4 


S-R3 


V G between positions 233 and 234 


Frame-sliift 


Same 


Truncation of CsaB 


G1024A 


csaB 5 


S-R4 


A A at position 203 


Frame-sliift 


Same 


Truncation of CsaB 


G1024A 


csaB 6 


S-R5 


C385T 


Nonsense 


ND 


Truncation of CsaB 


G1024A 


csaB 7 


S-R6 


V GC™ between position 219-220 


Frame-sliift 


Same 


Truncation of CsaB 


G1024A 


csaB 8 


S-R7 


T274G 


Y92D 


ND 


Noncliarged to cliarged (-) 


G1024A 


csaB 9 


S-R8 


A C at position 129 


Frame-sliift 


ND 


Truncation of CsaB 


WT 


csaB 10 


S-R9 


C254T 


T85M 


ND 


Polar liydropliilic tononpolar liydropliobic 


WT 


csaB 1 1 


JB220 parent 


WT 


None 


WT 




WT 


csaB 12 


J-Rl 


G at position 573 to T 


R191S 


Same 


Cliarged (+) to noncliarged 


WT 


csaB 13 


J-R2 


V G between positions 233 and 234 


Frame-sliift 


ND 


Truncation of CsaB 


WT 


csaB 14 


J-R3 


G421T 


Nonsense 


ND 


Truncation of CsaB 


WT 


csaB 15 


J-R4 


C 292T 


Nonsense 


ND 


Truncation of CsaB 


WT 


csaB 16 


J-R5 


G796C 


G266R 


ND 


Noncliarged to cliarged (+) 


WT 


csaB 17 


J-R6 


T258G 


S86R 


ND 


Noncliarged to cliarged (+) 


WT 


csaB 18 


J-R7 


V G between positions 233 and 234 


Frame-sliift 


ND 


Truncation of CsaB 


WT 


csaB 19 


J-R8 


C767A 


Nonsense 


ND 


Truncation of CsaB 


WT 



V and A indicate insertion and deletion of base(s) respectively at the indicated position from the start of the gene. ND- Not Done. 
WT- Wild type; G1024A indicates a SNP in BAS3946 at position 1024 from the start of gene from G to A. 



Bishop-Lilly et al. Virology Journal 2012, 9:246 
http://www.vi rologyj.conn/content/9/1 7246 



Page 7 of 14 



Table 4 Strains and plasmids used in this study 



Strain 


Description 


Source/Reference 


Escherichia coli SMIO 




[39] 


Bocillus onthfQcis 34F2 


bterne pxu i pxuz 


Hanna lab 


Bacillus anthracis BA663 


7702 Steme pXOl+pXOZ 


Koehler lab 


Bocillus onthrocis JB220 


7702 ABAS2245 


Calendar lab 


Bacillus anthracis S-Rl through S-R9 


Spontaneous AP50R mutants 


This study 


Bacillus anthracis J-Rl through J-R8 


Spontaneous AP50R mutants 


This study 


Bacillus anthracis BAP350 


BA663 AcsaB 


This study 


Bacillus anthracis BAP356 


BA663 

{/\rDU\i-D)-CSaD VLiZ34 


This study 


bacillus antnracis DArooo 


dAooo 

(AP50R-1 )-cs6/5 VG46 


This study 


Bacillus anthracis BAP41 1 


BA663 (AP50R-6)-csa8 VGCTO219 


This study 


Bacillus anthracis BAP435 


BA663 (AP50R-4)-csa8 AA203 


This study 


Bacillus anthracis BAP503 


BA663 csaB H270A 


This study 


Bacillus anthracis BAP533 


BAP350 csaB A693C creating an Xmal site 


This study 


Bacillus anthracis BAP553 


BA663 ABAS3946 


This study 


8cic/7/us anthracis BAP560 


BA663 with BAS 3946 G1024A mutation 


This study 


8cic/7/LiS anthracis BGBN 001 


BAP350/pBGBN1003 


This study 


Plasmid pRP1028 


Mutant construction 


This study 


Plasmid pSS4332 


Mutant construction 


This study 


Plasmid pSW4 


E. coli- B. anthracis shuttle vector 


[37] 


Plasmid pBGBNlOOl 


pSW4 NotI 


This study 


Plasmid pBGBN1002 


pSW4 oriT 


This study 


Plasmid pBGBN1003 


pSW4 oriT csaB 


This study 


Plasmid pG0v4 


pUC19 Amp^ Kan^ 


Gene Oracle 


Plasmid pNV36 


pG0v4 ABAS3946 


Gene Oracle 


Plasmid pNV37 


pG0v4 BAS3946 G1024A 


Gene Oracle 



V-insertion and A-deletion of base (s) at any given position from tine start of tlie gene; e.g., csaB VG46 indicates insertion of G base after position 46 (between 46 
and 47) from tlie start of the gene, and csaB AA203 means deletion of A at position 203 from the start of gene. 



other mutations that may be present in the original 
AP50^ mutants, we reconstructed four of the point 
mutations in a different genetic background, in Sterne 
strain 7702, which was not previously exposed to APSOc 
(see Methods). All four mutants (VG46, AA203, 
VGCTTA219, VG234) gave rise to mucoid colony 
morphology and phage resistant phenotype similar to 
the original spontaneous mutants, indicating that csaB 
indeed is involved in phage sensitivity. In addition to 
these mutants, we created a scarless, markerless in- 
frame deletion of csaB in 7702 background, which also 
showed mucoid and phage resistant phenotypes 
(Figure 3). 

Functional CsaB is needed for phage sensitivity 

Having shown that inactivation of csaB results in phage 
resistance, we wished to rule out the possibility that this 
result is due to a polar effect of the csaB mutations on 



downstream gene(s) or due to an alteration of regulatory 
elements in the deletion mutant. csaA and csaB form an 
operon upstream of sap and eae (Figure 2a), which have 
been shown to be involved in the formation of the cell sur- 
face associated S-layer (Sap) or an alternate S-layer (EAl) 
respectively [28]. In order to test the hypothesis that a 
functional CsaB is needed for phage sensitivity, we sought 
to inactivate the protein by mutating a critical amino acid 
residue while retaining the full-length protein and all the 
potential upstream and downstream control elements. A 
protein sequence alignment of 163 known orthologs 
revealed several conserved residues within CsaB protein 
(Figure 4 and Additional file 1: Figure SI) of which the 
histidine at position 270 is 100% conserved. We mutated 
histidine 270 to an alanine, and the resulting mutant was 
found to be mucoid and phage resistant (Figure 3), sug- 
gesting that a functional CsaB protein is necessary for 
phage infectivity, and inactivation of the CsaB protein 
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Figure 3 Phage sensitivity spot test on Sterne 7702 and various 
derivatives. For each row the left spot was spotted with medium 
(M), the middle spot received AP50c (A) and the right spot was 
spotted with y (G). The clearing in the bacterial spot indicates 
infection and killing of bacteria. In rows 3 through 6, the y spots 
completely cleared the bacterial spots. The various strains are as 
follows: 1) and 2) BA663 (wild type); 3) BAP350 {AcsaB); 4) BAP356 
{csoB VG234); 5) BAP366 {csoB VG46); 6) BAP41 1 {csaB VGCTO219); 
7) BAP435 icsaB AA203); 8) BAP503 {csaB H270A); 9) BAP533 {csaB 
A693C); 10) BAP553 (ABAS3946); 1 1) BAP560 (BAS3946 G1024A); 12) 
BGBN 001 (Aaa5/pBGBN1003). 



either by truncation or by mutating a critical amino acid 
residue gives rise to phage resistance. 

Complementation of csaB mutants with wild type copy of 
csaB 

We took two parallel approaches to test whether com- 
plementation of AcsaB mutant with a wild type copy of 
csaB restores phage sensitivity. First, we introduced a 
copy of csaB in its native locus in the csaB deletion 
strain. In this case, however, we engineered the gene to 




710 



715 



720 



position 



Figure 4 Amino acid sequence alignment of csaB orthologs. 

BAS0840 was aligned with the Pfam PF04230 seed alignment, and 
percent conservation was determined as described in Methods. Only 
the alignment of the amino acids in the vicinity of histidine 270 of 
CsaB protein is shown. The full-length alignment is provided in 
Additional file 1: Figure SI. 



carry a silent mutation (A693C) creating an Xmal site in 
order to distinguish this strain from the phage-sensitive 
parent strain. The resulting csaB-Xmal mutant strain 
was verified by PGR to carry the Xmal mutation 
(Figure 5a) and was found to be non-mucoid and phage 
sensitive. Thus, restoring the wild type allele in its native 
locus restored the phage sensitive phenotype. Second, 
complementation was done in trans by cloning the wild 
type csaB on a multicopy plasmid (Figure 5b) and 
transferring the recombinant plasmid (pNBGDlOOS) 
into the AcsaB strain via conjugation. In this plasmid, B, 
anthracis pagA promoter was used to drive the constitu- 
tive expression of the cloned csaB, The resulting recom- 
binant strain was found to be non-mucoid and phage 
sensitive (Figure 3). 

Transmission electron microscopic analyses of csaB 
mutants and complemented strains 

In an earlier study we found that the spontaneous AP50^ 
mutants produced an extracellular matrix and that AP50c 
phage failed to attach to the mutant bacterial surface 
[23]. In the current study, we confirmed these observa- 
tions by transmission electron microscopy (Figure 6a-d). 
As evident from these images, AP50c phage particles do 
not attach to the surface of the AcsaB mutant and conse- 
quently fail to infect. Phage adsorption and infectivity 
were regained when AcsaB was complemented in cis or 
in trans with the wild type copy of csaB (Figure 6e-h). 
The extracellular secretion is seen only in the mutants 
and is absent in the wild type and the complemented 
mutants. The TEM images also revealed that there may 
be a defective or improperly assembled S -layer, or the 
S -layer may have been masked by the presence of the 
extracellular matrix. Further confirmation of the adsorp- 
tion defect and reversion of this defect in complemen- 
ted strains was obtained in phage adsorption assays. In 
the csaB mutants (BAP350, BAP356, BAP366, BAP411, 
BAP435, BAP503) --62-100% of the phage particles were 
left in the supernatant whereas in wild type (BA663), 
AcsaB complemented strains (BAP533 and BGBNOOl) 
and the BAS3946 mutants (BAP553, BAP560) most of 
the phage (>99%) particles were adsorbed and very few 
particles (-0.1-0.65%) were left in the supernatant. 
Thus, the defect in phage adsorption was rectified in 
the complemented strains. 

csaB mutants are sensitive to y phage 

Earlier studies have shown that the cell surface protein 
GamR, which carries a LPXTG motif, is the cellular re- 
ceptor of the y phage [12]. We asked whether the cell 
surface-anchoring defect associated with csaB mutants 
affects GamR in any way, thereby preventing y adsorp- 
tion and subsequent infection. Consistent with earlier 
findings [30], all the mutants and the complemented 
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Figure 5 Complementation of csaB mutation, a) PGR verification of Xma\ mutant csaB. The wild-type sequence of csaB was restored in tlie 
AcsoB strain witli tine addition of a silent mutation creating an Xma\ site (see Methods). The region was amplified by PGR, and PGR products were 
purified and digested with Xmo\. The lane labeled BA663 indicates the csoB PGR product from wild type strain (7702) after digestion with Xmo\, 
whereas the lane labeled BAP533 is the PGR product from the mutant showing the two bands after digestion with Xma\. b) Genetic map of csaB 
complementing plasmid pNBGDlOOB. pNBGDlOOS is a shuttle vector derived from pSW4. The locations of the Gram positive and Gram negative 
replicons (from pUBllO and GolEl, respectively), antibiotic resistance genes oodD and bio, origin of conjugational transfer oriT, constitutively 
expressed B. anthracis pagA promoter (labeled PA promoter) and the csaB gene downstream of the promoter are indicated. 



strains allow infection of y phage, indicating that muta- 
tions in csaB do not affect gamR expression or bacterial 
surface localization of GamR protein (Figure 3). Thus, 
C5<3^5-mediated resistance is specific to AP50 phage and 
is not a general phage resistance mechanism. 



Conclusions 

Phage resistance mechanisms are an integral part of the 
evolutionary arms race between phage and bacteria that 
has been waged for millions of years. Generally two 
types of bacterial defense mechanisms are recognized: 




Figure 6 Transmission electron microscopic images of wt and AcsaB mutants and the complemented strains, a) BA663 (csaB+); b) BA663 
+AP50c; c) BAP350 (AcsaB); d) BAP350 (AcsaB)+AP50c; e) BAP533 (csaB+Xmal); f) BAP533+AP50c; g) BAP350 (AcsaB)+ pBGBN 1003 (csaB+); h) 
BAP350 (AcsaB)+ pBGBN 1003 (csaB+)+AP50c. The scale bar is in the bottom right corner of each figure. 
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active (restriction/modification systems, CRISPR/c<3^5, 
and abortive infections) and passive host adaptations [2]. 
In the latter category, at least two, if not three, types of 
host adaptations interfering with the phage life cycle can 
be hypothesized: 1) phage entry, 2) phage growth after 
entry, and 3) potentially lysogeny can be considered to 
be the first line of defense against phage infection since 
lysogens acquire resistance rather quickly and at very 
high percentage upon infection with a turbid plaque 
APSOt phage [23]. Mutants of phage that make clear pla- 
ques arise in the phage population; these mutants were 
found to carry two mutations within the phage genome 
— one in a promoter region and another in a gene of 
unknown function [23]. Although the clear plaque 
mutants do not seem to regain the infectivity of a lyso- 
gen in the case of AP50 (Sozhamannan, unpublished 
data), spontaneous resistant strains to APSOc phage arise 
at very low frequencies in bacterial populations. In 
mutants that are bacterial surface modified, phage resist- 
ance could be caused by masking of the bacterial recep- 
tor or by bacterial receptor mutations. Other mutations 
might affect subsequent steps after the phage DNA 
enters the cell such as DNA replication, particle assem- 
bly, or even lysis of the host. In an effort to understand 
the passive defense mechanisms, we isolated spontan- 
eous phage resistant bacterial mutants. These mutants 
secreted an extracellular matrix that appeared to protect 
the bacterium from phage infection. We hypothesized 
that this extracellular layer might mask the actual cellu- 
lar receptor, thereby preventing phage adsorption, since 
it is a gain of new function due to a mutation and is ab- 
sent in the wild type strain. 

Genetic mapping of a mutation entails linking the 
phenotype to a genetic locus. Conventional genetic 
methods use phage-mediated transduction, transform- 
ation or conjugation, by which the mutated gene is co- 
transferred with a known marker and mapped with 
reference to the known marker. Recombinant DNA and 
cloning techniques have been used to isolate and identify 
the gene responsible for a given phenotype. These con- 
ventional techniques are time-consuming and are not 
feasible in many organisms due to the lack of well- 
defined genetic systems. Recent developments in 
second-generation sequencing technologies have allowed 
us to obtain whole genome sequences of mutants and 
compare the sequences to that of the wild type strain to 
determine the causal variations. By taking a whole gen- 
ome sequencing approach, we have traced the causal 
variation of the phage resistant phenotype to csaB, There 
are two caveats to this approach. First, the sequencing 
technology used here (Roche/454) is highly error-prone, 
resulting in approximately 50-100 false-positive variant 
calls per bacterial genome sequence. Empirically, we 
have determined that by applying a filter based on the 



coverage and percentage of variant reads in a given pos- 
ition, positions with 75% or above tend to be true posi- 
tive variations. We have used a second filter in this 
study: a biological replicate strategy for ascertaining 
causal variations. A second caveat to this approach is 
that we assumed that the phenotype is due to a mono- 
genic trait and that there is a single pathway to give rise 
to that phenotype. However, despite these limitations, in 
the current study, we were able to correctly associate the 
causal variation of the phage resistance phenotype to 
csaB and to further validate this prediction by recon- 
structing various csaB mutants and showing the linkage 
between phage resistance and csaB mutations. 

Mesnage et al. have characterized the csaAB operon of 
B. anthracis extensively [28]. These authors have pro- 
posed a role for csaB in cell surface anchoring of S -layer 
homology (SLH)-containing proteins via pyruvylation of 
cell wall polysaccharides. Consequently, mutants lacking 
csaB appear to be defective in cell surface anchoring of 
various SLH proteins including EAl and Sap. In support 
of this idea, in csaB mutants. Sap and EAl are conspicu- 
ously absent from the cellular fraction and are primarily 
found in the secretome. These authors also observed 
additional aberrant cellular morphologies that were 
replicated in this study, including the accumulation of 
cell wall material on the cell surface. These authors have 
attributed this phenotype to the lack of autolysis func- 
tion in csaB mutants [28]. Recently, the role of yet an- 
other gene, tagO, has been implicated in S -layer protein 
assembly via binding to secondary cell wall polysacchar- 
ide [31]. Given the known function of CsaB in the pyru- 
vylation pathway of polysaccharides, its most probable 
role in AP50 life cycle may be indirect: in the proper cell 
surface anchoring of phage receptor protein. The most 
probable candidates for the bacterial receptor may be 
Sap and EAl. It is possible that Sap or EAl is the 
AP50 receptor, and csaB mutants fail to support AP50 
adsorption because of the defect in proper cell surface 
anchoring of the S -layer. The transmission electron 
micrographs support the idea that the csaB mutants 
lack a well-organized S -layer. 

The spontaneous phage resistant mutants isolated in 
this study carry different types of mutations: SNPs and 
small insertions leading to nonsense codons, frame- 
shifts and non-synonymous changes. While a vast major- 
ity of these mutations cause CsaB truncation, some 
other mutants carry non-synonymous changes in the 
full-length protein. The essentiality of these amino acid 
residues in the overall structure and function of CsaB is 
not readily apparent since these mutations fall in regions 
of very low conservation in the protein alignment 
(Figure 4 and Additional file 1: Figure SI), but it is pos- 
sible that the individual amino acid changes may affect 
protein folding and therefore function. The AP50^ 
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mutants derived from 34F2 also carried a mutation in 
BAS3946 besides csaB, We have clearly established that 
the mutation in BAS3946 does not contribute to phage 
resistance in any way, and it is probably a 'hitch-hiker' 
mutation that was carried in the immediate parent prior 
to the selection for phage resistance. This result also 
highlights the power of the next generation sequencers 
for unbiased mapping of all variations arising in a given 
bacterial strain under selection pressure or mutagenesis, 
unlike conventional genetic techniques, which scan a fo- 
cused area of the genome for mutations. 

Another intriguing question that remains is why spon- 
taneous AP50 resistant mutants were not found in the 
actual phage receptor gene itself, assuming Sap or EAl 
is the AP50 phage receptor. It is possible that there are 
differences in the mutability or essentiality of the genes, 
although we have already shown that sap is not an es- 
sential gene for B. anthracis viability (Plant et al, unpub- 
lished data). Further characterization of additional genes 
involved in APSOc life cycle identified by this study is 
underway. The fact that y phage and AP50 use different 
receptors is useful in phage therapy experiments, since 
development of resistance against both phages would be 
kept under check when used in combination. 

Methods 

Bacterial strains, plasmids and oligonucleotide primers 

Strains and plasmids used in this study are listed in 
Table 4. The three Sterne strains used in this study were 
34F2 (obtained from Philip C. Hanna), 7702 (obtained 
from Theresa M. Koehler) and JB220 (obtained from 
Richard Calendar). Strain 7702 is designated as BA663 
in this study. Derivatives of 34F2 that were spontan- 
eously resistant to AP50 are designated S-Rl through 
S-R9, and those of JB220 are designated J-Rl through 
J-R8. When necessary, strains were grown in the pres- 
ence of spectinomycin (100 (ig/ml for E, coli; 250 (ig/ml 
for B. anthracis) or kanamycin (20 (ig/ml for E, coli; 25 
(ig/ml for B, anthracis). The primers used in various 
plasmid constructions and Sanger sequencing of PGR 
amplicons are listed in Additional file 2: Table SI. 

Genome sequencing and variation detection 

Whole-genome sequencing (WGS) was performed in the 
Genome Sequencer FLX (454 Life Sciences/Roche) using 
FLX or FLX Titanium reagents according to the manu- 
facturer's protocols and instructions [32]. Signal proces- 
sing of FLX data was performed on the sequencer itself, 
whereas signal processing of Titanium data was per- 
formed off-rig on a Linux cluster of 10 nodes connected 
via gigabit ethernet. Each node contained eight 64-bit 
processing cores running at 2.3 GHz with 8 GB of RAM. 
De novo or reference-guided assembly of either FLX or 
Titanium sequences, preliminary annotation of the 



genome using the DIYA pipeline [33] and identification 
of variations were performed using the above-described 
computer cluster and Sterne (NC 005945) as reference 
genome as described [27]. A variant flagged by GS Refer- 
ence Mapper was defined as high confidence if it met 
the following criteria: a) at least three non-duplicated 
reads showed the same variation; b) at least one read in 
each of the forward and reverse orientations showed this 
variation, unless there are at least five reads with quality 
scores over 20 (or 30 if the difference involves a 5-mer 
or higher). A more detailed description of differences is 
available in the Genome Sequencer Data Analysis Soft- 
ware manual. In this study, a high quality variation is 
defined as a high confidence (as defined by 454/Roche) 
and high concordance (>75% reads have the variation) 
variation. 

Culture methods and phage assays 

All bacterial strains were grown in Luria Broth or Brain 
Heart Infusion medium, and the phage assays were per- 
formed according to published protocols. Phage spot tests 
were performed as follows: first, triplicate spots of 50 \A 
of an overnight culture were made on large phage assay 
agar plates and allowed to dry. Five (il of phages (-10^ 
PFU) AP50c and y was spotted on two bacterial spots 
while the third spot received 5 (il of medium. The spots 
were dried at room temperature and the plates were sub- 
sequently incubated at 37°C overnight. Phage adsorption 
assay was performed as follows: One ml of log phase wild 
type and mutant B, anthracis cells was spun and the pel- 
let was resuspended in 100 \A of AP50c phage and incu- 
bated for 30 minutes at room temperature. After addition 
of 900 (il of phage assay broth, the bacteria-phage com- 
plex was filtered through a 0.45 (im syringe filter. The 
plaque forming units (PFU) in the filtrate was titred using 
7702 as the indicator bacterium. The percentage adsorp- 
tion was calculated against the same volume of unad- 
sorbed phage filtered through the syringe filter. 

Construction of mutants 

Two plasmids were engineered to facilitate allelic ex- 
change in B, anthracis, as improvements to a published 
system [29]. Details of the improvements will be pub- 
lished elsewhere. Briefly, pRP1028 and pSS4332 serve 
the functions of pBKJ236 and pBKJ223, respectively, and 
contain many of the same features. In pRP1028, the 
erythromycin resistance of pBKJ236 was replaced with 
spectinomycin resistance, and in pSS4332, the tetracyc- 
line resistance of pBKJ223 was replaced with kanamycin 
resistance. Allelic exchange constructs derived from 
pRP1028 were transferred to B. anthracis, and plasmid 
integrants were isolated following a temperature shift 
[29]. Introduction of pSS4332 into the integrant led to 
I-<SceI-mediated cleavage of the integrated plasmid. 
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stimulating the second crossover event. The presence of 
the desired changes was confirmed by PGR and 
sequencing. 

Specifically, to reconstruct the point mutations found in 
the spontaneous AP50^ mutants, the mutations and ap- 
proximately 500 bp of flanking homology on each side 
were PCR-amplified from the B. anthracis strains and 
cloned into pRP1028. To delete csaB, standard cloning 
methods were used to engineer a construct with the se- 
quence GTG CGG CCCGGG GGA TCT TAA, represent- 
ing START+one co&on^-Xmal site+two codons+STOP, 
along with approximately 500 bp of flanking homology. 

To restore the wild-type allele to the csaB deletion 
strain, a derivative of pRP1028 was engineered with the 
csaB sequence with the adenosine at position 693 chan- 
ged to cytosine (creating the Xmal site CCCGGG without 
changing the amino acid sequence), along with approxi- 
mately Ikb of flanking homology. The change was con- 
firmed by PGR and sequencing, as well as by digestion of 
the PGR product with Xmal (Figure 5a). 

To identify conserved residues of CsaB, the nucleotide 
sequence (BAS0840) was aligned with the 163 proteins in 
the seed alignment of Pfam family PF04230 (polysacchar- 
ide pyruvyl transferase) [34] using ClustalW [35]. The 
resulting alignment was processed with WebLogo [36] to 
calculate the consensus sequence and percent conserva- 
tion at each position (Figure 4 and Additional file 1: 
Figure SI). The histidine at position 270 of BAS0840 was 
found to be 100% conserved in these proteins. Standard 
cloning techniques were used to engineer a derivative of 
pRP1028 in which the histidine codon (CAT) was chan- 
ged to alanine (GCT) (Figure 2), flanked by approximately 
850 bp of homology. The presence of the H270A change 
was confirmed by PGR and sequencing. 

DNA sequences used to generate ABAS3946 and 
BAS3946 G1024A were synthesized (Gene Oracle, 
Mountain View, CA). For ABAS3946, the sequence 
included the first 6 codons of the gene, followed by the 
last 5 codons of the gene and the STOP codon, all 
flanked by approximately 500 bp of homology. For 
BAS3946 G1024A, the sequence included the mutation 
flanked by 500 bp of homology. In each case, DNA was 
cloned into pRP1028, allelic exchange was performed as 
described above, and the presence of the desired changes 
was confirmed by PGR and sequencing. 

Construction of pNBGD1003 and complementation of AcsaB 

The starting plasmid for pNBGD1003 was pSW4 [37] 
which was modified by inverse PGR to include a Notl- 
Ndel site on either side of the pagA promoter using pri- 
mers 1004 and 1005. The pagA promoter was amplified 
using primers 1006 and 1026, and fragments were 
digested with Notl and Ndel and cloned, giving rise to 
pNBGDlOOl. Plasmid pNBGD1002 was derived from 



pNBGDlOOl by introducing an origin of transfer, oriT, 
as a Notl fragment at the Notl site. The oriT fragment 
was PGR amplified from pIN297, which is a derivative of 
pPL2 [38] using primers 1064 and 1065. Plasmid 
pNBGD1003 was derived from pNBGD1002 by clon- 
ing a Pmll' Bglll fragment carrying the csaB gene 
sequences into Ndel digested, blunted and BamHl 
digested pNBGD1002. Plasmid pNBGD1003 was elec- 
troporated into SMIO and subsequently transferred 
into B. anthracis AcsaB mutant by conjugation, 
selecting for polymyxinB and kanamycin resistance as 
described previously [29]. The resulting ex-conjugants 
were purified and checked for phage adsorption and 
infection. 

PGR and Sanger verification of the csaB and BAS3946 
genes from APSO'^ mutants 

The DNA sequences spanning the csaB (1385 bp) and 
BAS3946 (1445 bp) genes were PGR amplified from wild 
type and various mutants (both 34F2 and JB220 deriva- 
tives), and the PGR amplicons were purified and 
sequenced on both strands using primers indicated in 
Table S-1 (BGBN14 through BGBN21). The resulting 
nucleotide sequences were aligned to the wild type 
sequences of these two fragments to determine the mu- 
tational changes. 

Transmission electron microscopy 

Bacterial pellets were fixed in 2.0% glutaraldehyde in 0.1 
M sodium cacodylate pH 7.4 overnight at 4°C. Following 
buffer rinse, samples were post- fixed with 1% osmium 
tetroxide in 0.1 M sodium cacodylate for 1 hr on ice in 
the dark. After a brief distilled H2O rinse samples were 
placed in 2% aqueous uranyl acetate (0.22 um filtered) 
for 1 hour at room temperature in the dark. Following 
en-bloc staining, cells were dehydrated through a graded 
series of ethanol to 100% ethanol, transferred through 
propylene oxide, embedded in Spurr s low viscosity resin 
(Polysciences), and cured at 60°C for two days. Sections 
were cut on a Riechert Ultracut E with a Diatome Dia- 
mond knife. Eighty nm sections were picked up on for- 
mvar coated 1x2 mm copper slot grids and stained 
with uranyl acetate followed by lead citrate. Grids were 
viewed on a Hitachi 7600 TEM operating at 80 kV and 
digital images captured with an AMT 1 K x 1 K CCD 
camera. 

Additional files 



Additional file 1: Figure SI. Sequence alignment of csaB and 
orthologs. BAS0840 was aligned with the Pfam PF04230 seed 
alignment and percent conservation was determined as described 
in Methods. 
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Additional file 2: Table SI. List of primers used ir^ this study. 
Restriction sites are ur^derlined. 
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